Lost at Sea 



Steven R. Finch and John A. Shonder 

November 23, 2004 

Abstract. What is the path of minimum expected length for escaping 
a planar convex region f2? We rigorously obtain best 2-segment and 3-segment 
solutions when is an infinite strip, and numerically examine 2-segment solu- 
tions when f2 is a disk. 

A swimmer is lost in a dense fog at sea. She knows that the sea is a planar infinite 
strip 

{(x,y) E R 2 : < x < 1} 

of unit width. Assume that the x-coordinate of her initial position is uniformly dis- 
tributed on the interval [0, 1]; the distribution of the y-coordinate can be arbitrary. 
Assume as well that her initial orientation (the angle between her initial velocity vec- 
tor and the x-axis) is uniformly distributed on [— n, n] and that her speed is constant. 
What escape trajectory should the swimmer follow that minimizes her expected time 
to reach either shore? 

The min-max analog of this problem was solved long ago (see PP for a survey). 
Zalgaller 00] proposed a heuristic solution of the above min-mean problem. Shonder 
|1] independently determined the precise 2-segment solution of least expected escape 
time. (We agree that a fc-segment path is a continuous, piecewise linear curve con- 
sisting of < k pieces.) The details underlying Shonder's computation appear here for 
the first time. We find a 3-segment path that improves slightly upon the 2-segment 
path and also demonstrate that Zalgaller's solution is far from optimal. 

Next, we examine a different sea: the disk 

{{x,y) G M 2 : x 2 + y 2 < l} 

of unit radius. Assume that the starting point is uniformly distributed on the disk 
and that we seek (as before) to minimize the expected time to reach the boundary. 
Intuition suggests that the swimmer should follow a 1-segment path. If this is true, 
then the mean length of an arbitrary escape trajectory is > 8/(37r) = 0.8488263631... 
A calculus-of- variations proof of this general inequality is not known. We can 
confirm this only numerically for 2-segment paths. 

In the final section, there appears a first attempt at proving the above inequality 
under special circumstances. The technique is due to Gevirtz jH] and yields results 
when the escape path curvature is sufficiently small. 

°Copyright © 2004 by Steven R. Finch. All rights reserved. 
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Figure 1: Five cases for the infinite strip. 

0.1. Infinite Strip: 2-Segment Scenario. Without loss of generality, let the 
initial position be (x, 0) and the initial orientation be 9, where < x < 1 and 
< 9 < 7i. Fix a distance r > and an angle < a < tt. As shown in Figure 1, 
there are five distinct cases to consider. 

Case 1: the swimmer reaches the right-hand shore before traveling a distance 
r, which happens when x + r > 1 and 9 < arccos((l — x)/r). The path length is 
(1 -x)/cos(0). 

Case 2: the swimmer travels a distance r, does not reach the right-hand shore, 
pivots an angle a, and then travels until the left-hand shore is reached. This happens 
when 

{[x + r > 1 and arccos((l — x)/r) < 9] or [x + r < 1]} and {9 < | — a} . 

The path length is r + (x + r cos(0)) / cos(# + a). 

Case 3: the swimmer travels a distance r, does not reach a shore, pivots an angle 
a, and then travels until the right-hand shore is reached. This happens when 

{[x + r > 1 and x — r < and arccos((l — x)/r) < 9 < n — arccos(a;/r)] 
or [x + r < 1 and x — r < and 9 < n — arccos(x/r)] 
or [x + r > 1 and x — r > and arccos((l — x)/r) < 9] 
or [x + r < 1 and x - r > 0]} and {f-a<6><^-a}. 

The path length is r + (x — 1 + r cos(0)) / cos(# + a). 
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Case 4: the swimmer travels a distance r, does not reach the left-hand shore, 
pivots an angle a, and then travels until the left-hand shore is reached. This happens 
when 

{[x — r < and 9 < n — arccos(x/r)] or [x — r > 0]} and {4f — a < 9} . 

The path length is r + (x + r cos(6>)) / cos(# + a). 

Case 5: the swimmer reaches the left-hand shore before traveling a distance r, 
which happens when x—r < and ir — arccos(x/r) < 9. The path length is —xj cos(0). 

We've restricted attention to < 9 < tt for simplicity's sake. The set of es- 
cape paths corresponding to — 7r < 9 < is obtained from the set of escape paths 
corresponding to < 9 < n via reflection across the horizontal axis, followed by 
reflection across the vertical axis. This composite transformation can be written as 
(x,9) i-> (l-x,9-n). 

If it is assumed that r > 1 and < | — a < arccos (1/r), then only Cases 1, 3 and 
5 enter into the calculations of the expected path length: 



i 

1 

7T 



7r— arccos 



1 - x Jn f ( x - 1 + r cos(#) \ Jn f x 

— d9 + / r + )d0- / --rd9 

cos(0) J V cos(6» + a) J J cos(6>) 

arccos (^-^) T— arccos 



dx. 



All of these integrals can be evaluated in closed form. Hence the expected path 
length, multiplied by ir, becomes 



1 + cos a) ^— + \J r 2 — 1 — r j — (2 + cos a) arccos ^-^ + - 

i ( r~o 7 \ 1 21 /-, \ • 2 -i / Vr 2 — 1 sin a — cos a; \ 

m yr z -1 + rH — r m(l — cos a) sin a + r sin am + 

V / 2 \ rsina / 

- (r 2 sin 2 a + l) In ^V?" 2 — 1 cos a + sin a + r j + 

- (r 2 sin 2 a — l) In (^Jr 2 — 1 cos a — sin a + r j — 

- (r 2 sin 2 a + l) In —y/r 2 — 1 cos 2 a — ^sin a + V?" 2 — 1 — r^j cos a; — sin a + r J — 

- (r 2 sin 2 a — l) In —y/r 2 — 1 cos 2 a + ^sina — Vr 2 — 1 + cos a + sin a + rj . 

Minimizing this expression gives r = 1.0432668686... and a = 1.3734935859... ~ 78.7° 
as the optimal parameter values for the 2-segment scenario. Therefore the least 
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Figure 2: Four sample realizations with r = 1.043, a = 78.7°, x = 0.1,0.3,0.5,0.9 
and 9 = 80°, 160°, 50°, 140° (involving three Cases). 
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Figure 3: Two subcases of Case 3. 

expected path length is 0.8869669056.... See Figure 2 for several sample realizations 
(with fixed r and a). 

Here is a new problem. For each (r, a) , there exists a probability distribution of 
escape times obtained by sampling (x, 9) uniformly from [0, 1] x [0, n]. Compute the 
median of each distribution and then determine (r, a) corresponding to the least me- 
dian. Monte Carlo simulation suggests that a solution is r = oo (a can be arbitrary), 
consistent with the idea that medians are less sensitive to outliers than means. Define 
a random variable q to be (1 — x)/ cos(6 l ) if 9 < rr/2 and —xj cos(9) if 9 > ir/2. The 
path length is the median of q, which is estimated to be 0.78.... 

0.2. Infinite Strip: 3-Segment Scenario. We assume that r > 1 and < 
| — a < arccos(l/r) as at the end of the preceding section. Cases 1 and 5 thus 
remain unchanged and only Case 3 needs to be refined. Fix a second distance s > 
and a second angle < (3 < ir. See Figure 3. 

Subcase 3.1: the swimmer travels a distance r, does not reach a shore, pivots 
an angle a, and then reaches the right-hand shore before traveling a distance s. This 
happens when 

arccos((l — x)/r) < 9 < 9 (x, r, s, a) 
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where #0 is a solution of the equation 

1 — x — r cos(#) + s cos(# + a) = 0; 
an exact formula for Q Q is 

1 -x \ 

^Jr 2 + s 2 — 2 r s cos(a) J 

The path length is r + (rr — 1 + r cos(0)) / cos(# + a). 

Subcase 3.2: the swimmer travels a distance r, does not reach a shore, pivots 
an angle a, travels a distance s, does not reach a shore, pivots an angle /3, and then 
travels until the right-hand shore is reached. This happens when 

(x, r, s, a) < 9 < n — arccos(:r/r) 

and (3 satisfies other constraints (for example, if f3 is too small, then the 3-segment 
might cross itself and reach the left-hand shore rather than the right-hand shore). 
We will not attempt to prescribe these additional inequalities. The path length is 
r + s + (1 — x + s cos(6 l + a) — r cos(0)) / cos(# + a + (3). 



. , ssin a 

B Q = arcsm | —^^^^^=^^^^^= ) + arccos 
yr 2 + s 2 — 2 r s cos(o;) 



Under the proper circumstances, therefore, the expected path length is 
1 — x ,„ f ( x — 1 + r cos(#) 



cos(fl) 



d9 + 



r + 



7T— arccos 



/ ( 

80 



r + s + 



arccos ( ) 

1 — x + s cos(6> + a) — r cos(#) 
cos(# + a + (3) 



cos(6> + a) 



d9+ 



TT 

)«,- J 



X 



tt— arccos 



3) 



cos(#) 



d9 



dx. 



A closed-form expression is unlikely here because of the presence of 9 in the middle 
two integrals: only integration with respect to 9 seems feasible. Let 

k — \Jr 2 + s 2 — 2 r s cos(a) , p = arcsin (s sin(a)//c) , 



u\ — r sin(o;), u 2 = s sin(/3) — r sin(o; + 
t»i = r(l + cos (a)), f 2 = r + s + scos(/3) — rcos(o; + 
£j = a + arccos ((1 — x)/r) , £ 2 = a + + P + arccos ((1 — x)/k) , 
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7] l = a + p + arccos ((1 — x)/k) , r] 2 = a + /3 + 7T — arccos(x/r) 
and, for j = 1,2, define 



Pi = U 3 ln 



cos(^) 


(-iY(l-x) 


cos(^) 


2 



sin(^) (sin(^) - l) + sin(^) - 1 j J 3 3 ' 



Then the expected path length, multiplied by ir, becomes 

i 

ln (^\/r 2 — 1 + r j + r — \/r 2 — + J (pi + P2) dx. 



This can be minimized by use of Leibniz's rule for differentiation under the inte- 
gral sign. We obtain r = 1.0255050653..., a = 1.4909825316... w 85.4°, s = 
0.5306340577... and (3 = 2.7495709960... w 157.5° as the optimal parameter values 
for the 3-segment scenario, with least expected path length 0.8835534788.... 

A third subcase, for which 7r < (3 < 2n, awaits examination. Either shore could 
be the final resting place. 

0.3. Infinite Strip: Zalgaller's Proposed Solution. Consider the path in Fig- 
ure 4, where the points A = (0,0), B = (0.814,0), C = (0.8460,0.0005), D = 
(1.3017,0.0151), E = (0.814,1) are joined by line segments, with the exception of 
points B, C which are joined by a tiny circular arc of radius 1, center E and angle 
0.032. Zalgaller j2J E] claimed that, if the swimmer follows this path, then her ex- 
pected time to reach either shore is approximately minimal. We examined his claim 
by computing best 2-segment and 3-segment fits to Zalgaller's path, and then calcu- 
lating the expected escape time via our formulas. For example, in the 2-segment fit, 
we obtained r = 1.3017, a = 64.3° and hence the expected escape time is 0.9188. 
This is somewhat consistent with Zalgaller's estimate 0.9523. More importantly, how- 
ever, Shonder's optimal 2-segment path jl] has expected escape time 0.8870 which is 
considerably less than 0.9523. There must be an error somewhere in the details of [2]. 
We aim someday to better understand Zalgaller's elaborate construction, in the hope 
that his procedure (once corrected) will lead to escape trajectories that outperform 
even our optimal 3-segment path. 

0.4. Disk: 2-Segment Scenario. Without loss of generality, let the starting 
point be (x, 0) and the initial orientation be 8, where < x < 1 and —71 < 9 < 71. 
Fix a distance < r < 2 and an angle < a < 71. As shown in Figure 5, there are 
two distinct cases to consider. For convenience, define 

1 — x — r \ , . / r ' 



(p(x, r) = arccos ( j , r) = arccos ^ j 
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Figure 4: Zalgaller's escape path. 



Figure 5: Two cases for the disk. 
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Case 1': the swimmer reaches the shore before traveling a distance r, which 
happens when x > |r — 1| and — tp < 9 < ip, or when x < r — 1. The path length is 
q = — xcos(9) + \/l — x 2 sin 2 (0). 

Case 2': the swimmer travels a distance r, does not reach the shore, pivots an 
angle a, and then travels until the shore is reached. The path length is r + s, where 



y cos(a + uo) + y 1 



y 2 sin 2 (a 



y = \/x 2 + r 2 + 2xr cos(#), 



arcsm 



a; 



x sin(#) 



7r — arcsm 



x sin(9) 



if (x > r and — ip < 9 < ip) or x < r, 
if x > r and (9 < —ip or ip < 9) . 



The proof of Case 1' makes use of the Law of Cosines and the fact that (p exists if and 
only if x > \r — 1|. The Law of Sines appears in the proof of Case 2' when defining 
the angle u of an auxiliary triangle with side lengths r, x and y; the Law of Cosines 
is then applied to the triangle with side lengths 1, s and y to obtain s. 
The expected path length is therefore I + J, where 



7=1 

7T 



/ (x, r, a) xdx, 



J 



r-l 



— / g (x) xdx 



|r-l| 



7T 



-r+1 



if r > 1, 



1 

7T 



/i (x, r, a)idx if r < 1 



and 



/ = y q d9 + J{r + s)d9 + J{r + s)d6, g = J qd9, h = J{r + s)d9. 



Closed-form expressions are not possible here. Numerical computations confirm that 
8/(37r) is the value of / + J when a = ir, for any r > 0, and that otherwise I + J is 
strictly larger than 8/(37r). See Figure 6 for several sample realizations (with fixed r 
but variable a). 

We mention a corresponding min-median result: Monte Carlo simulation suggests 
that a solution is r = 2 (a can be arbitrary). The path length is the median of the 
random variable —xcos(9) + a/1 — x 2 sin(6') 2 , which is estimated to be 0.94.... 

These kinds of "lost at sea" problems deserve to be more thoroughly studied! The 
only other published reference we know regarding the min-mean problem is jjj- 
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Figure 6: Sample realizations with r = 0.5, x = 0.1,0.3,0.5, a = 110°, 160°, 20° and 
6 = 200 o ,80°,45°(all from Case 2'). 

0.5. Disk: Gevirtz's Technique. Let 7 : [0, 00) -^Cbea differentiable curve 
in the complex plane satisfying 7(0) = 0. Parametrize 7 by arclength s and let 
7'(s) = exp(i (f>(s)) , where we assume that <f>(s) = arg(7'(s)) is real, 0(0) = and 
|0( s )| < 0o < f • The function s 1— > A(s) = \"f(s) | is one-to-one and hence possesses 
an inverse A 1— > s(A). Let D denote the unit disk in C and define C(z) to be the circle 
of radius 1 centered at z; as an example, C(0) is the boundary of D. 
We wish to minimize the average 



where cr(z, 7) denotes the first (and only) value of s for which \j(s) — z\ — 1. (By 
symmetry of -D, we have substituted —2 for z.) Now, 




(7(2;, 7) dx dy 



D 



a(z, 7) dx dy 




D 







where fj,(s) is the area of the portion of D not enclosed by (7(7(5)). This identity 
follows via a Riemann sum argument. 
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Consider the family of circular arcs F = {.DnC(7(s)):0<s<oo}. If O is 
small enough, then no two distinct arcs in F can intersect. As a consequence of the 
above identity, if the curvature 4>'(s) is suitably small, then 



where s* satisfies |7(s*)| = 2. Note that s* depends on 7, which complicates any 
variational approach to this problem. Because A'(s) = cos (<j>(s) — arg(7(s))) and 
s(A) > A, we have 



as was to be shown. It is unclear what to do in the general situation for which distinct 
arcs in F are not necessarily disjoint. We also wonder if a similar technique exists for 
the infinite strip sea. 

0.6. Acknowledgements. We are grateful to Julian Gevirtz and Patrice Le Conte 
for their help. 
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